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The Ga vacancy mediated microstructure evolution of (Ga,Mn) As during growth and post-growth 
annealing is studied using a multi-scale approach. The migration barriers for the Ga vacancies and 
substitutional Mn together with their interactions are calculated from first principles, and temporal 
evolution at temperatures ranging from 200 to 350° C is studied using Lattice Kinetic Monte Carlo 
simulations. We show that at the typical growth and annealing temperatures (i) gallium vacancies 
provide the diffusion mechanism for substitutional Mn and (ii) in 10-20 h the vacancy mediated 
diffusion of Mn promotes the formation of substitutional clusters. Clustering reduces the Curie 
temperature (Tc), and therefore the Mn clustering combined with the fast interstitial Mn diffusion 
explains the experimentally observed twofold annealing behavior of Tc- 
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Understanding of microstructure evolution during 
growth and post-growth annealing is one of the key is- 
sues in materials science. The microstructure and its in- 
homogeneities largely determine the material properties, 
including the basic phase transition points for materi- 
als ranging from high temperature superconductors to 
diluted magnetic semiconductors P, • In the ongoing 
quest for room temperature semiconductor spintronics 
materials (Ga,Mn)As has been one of the main candi- 
dates since the observation of the relatively high Curie 
temperature Tc of 110 K Q. Typically the (Ga,Mn)As 
thin films, where Mn substitutionally replaces Ga atoms, 
are grown by means of low-temperature molecular beam 
epitaxy, which also leads to the formation of As antisites 
(AsGa) d and interstitial Mn (Mni) H 0. The sub- 
stitutional Mn (MuGa) act as acceptors, providing spin- 
polarized holes that mediate the ferromagnetic coupling, 
while the Asca and Mn; as donors hamper the ferromag- 
netism by compensating holes. Post-growth annealing of 
a few hours is an efficient method to remove the intersti- 
tial Mui and thus increase Tq However, ex- 
tended annealing at temperatures around 250°C for ten 
hours or longer reduces Tc again 0, 13 ' ^"^^ annealing 
at higher temperatures leads to a more swift lowering 
of Tc m. This twofold behavior clearly indicates that 
besides the Mui out-diffusion another yet unknown mi- 
crostructure evolution process takes place. At present, 
the nature of this process and the mechanism underlying 
this evolution is not understood. 

In this Letter we show that the mechanism behind the 
long-term microstructure evolution is gallium vacancy 
(Vca) mediated MuQa diffusion on the Ga sublattice. 
We also show that this diffusion leads to Mn cluster- 
ing, which reduces Tc [l(l ITU . Although the forma- 
tion of MuGa clusters has been shown to be energetically 



favorable [ill IT^l IT3 . ITHI |. it requires an abundance of 
mobile gallium vacancies. The recent discovery of rather 
high gallium vacancy (Vpa) concentrations in (Ga,Mn) As 
up to 10^^ cm"'' [l^ gives us good reason to consider 
this mechanism plausible, but the mobilities of Vca and 
MuQa in (Ga,Mn)As are unknown. We calculate their 
migration and binding energies from first principles and 
use them as input data in lattice kinetic Monte Carlo 
simulations of microstructure evolution. 
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FIG. 1: Mn/Ga migration path. The intermediate stages, 
where the Mn/Ga atom is in a planar configuration with sur- 
rounding As atoms 2, 4, and 5, and where the Mn/Ga atom 
passes through the (IlO) plane close to the Tca interstitial site 
are also given, and bonds in these configurations are denoted 
with the thick dashed lines. 

The vacancy mediated substitutional Mn migration 
over the Ga-sublattice consists of the following three 
steps: (1) the Mnca atom and a vacancy form a pair; 
(2) the MuQa and the vacancy switch places; and (3) the 
pair dissociates. This mechanism is henceforth called the 
pair formation-dissociation mechanism. The considera- 
tion of Mn diffusion via Vca requires the knowledge of at 
least the following three energy quantities: the migration 
energy for Vca, the migration barrier for the Mnca^Vca 
transition, and the MnGa~VGa interaction potential (i.e. 
binding energy). These energies are calculated using first 



principles methods. 

Total energy density-functional calculations of 
(Ga,Mn)As are performed using the supercell approach 
and the all-electron projector augmented-wave method 
together with the generalized gradient approximation 
(GGA-PW91) for exchange-correlation implemented in 
the VASP code [l^. The binding energies are calculated 
as the difference in total energy of the pairs located at 
nearest-neighbor and fifth-nearest-neighbor separations 
on the Ga fee sublattice in 64 atom supercells that are 
fully relaxed, as described in Ref. 0|. In the 107 atom 
supercell we use a 4 x 4 x 3 fc-point sampling mesh giving 
a similar fc-point density as that used in Ref. . 

As the result we find a binding energy of £^6 = 0.1 eV 
for the Mnoa dimer. We find further that the binding of 
a MuGa atom to a Mnca cluster can to good accuracy 
be estimated as n ■ Eh, where n is the number of nearest 
neighbor Mnca bonds formed. Consequently, the lowest 
energy microstructure can be expected to be composed 
of Mn clusters, i.e. a collection of Mn atoms occupy- 
ing nearest neighbor Ga sites, instead of separated single 
Mnca, or even of second phase precipitates for very high 
Mn concentrations. However, to approach the lowest en- 
ergy microstructure within reasonable time requires an 
efficient diffusion mechanism for the substitutional Mnca 
atoms. Therefore we consider the plausible Voa mecha- 
nism. However, in contrast to the MnGa~MnGa binding, 
we find that bringing a MnGa and VGa to neighboring 
sites is energetically unfavorable and costs 0.2 eV, i.e 
Eb = —0.2 eV. The negative binding energy means that 
the pair dissociation barrier will be 0.2 eV lower than 
the pair formation barrier. Nonetheless, such pairs may 
be formed kinetically if the formation barrier is reason- 
ably low. Thus the efficiency of the vacancy mediated 
diffusion is determined by the relation of the activation 
barriers for direct MnGa ^ VGa exchange and MuGa^ 
VGa pair formation-dissociation processes. 
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FIG. 2: Migration energy barriers for Ga (a) and Mn (b) 
difTusion vs. the lowest energy path (see Fig. ^ projection 
on the [110] axis. The crosses correspond to lowest energy 
points from the calculation with no relaxation, while the cir- 
cles indicate the saddle point energies from the calculation 
with constrained relaxation. The lines are only to guide the 
eye. 

The migration barriers are calculated using the drag 
method in a supercell of 107 atoms and one vacancy. 
We consider migration of Mn/Ga in a Mn rich environ- 
ment, i.e. migration in metallic (Ga,Mn)As. Thus, our 



supercells contain a Ga-VGa or a MnGa~VGa pair and 
one inactive substitutional MnGa at the largest distance 
from this pair. The Ga — > Vca and Mnca VGa jumps 
take place along a path in the (110) plane (Fig.^. First 
the migrating atom (Ga or MnGa) is dragged towards the 
VGa, and the lowest energy path is determined by calcu- 
lating several points on a grid around the assumed path 
(Fig. ^1 without relaxations. We assume that the jump 
process is symmetric, so we only need to calculate half of 
the path. The obtained migration path passes close to the 
Tca interstitial site, as shown in Fig. ^ During this pro- 
cess the migrating atom first breaks the bond to As atom 
1, followed by the breaking of bonds to atoms 4 and 5, 
thus only maintaining the bond to atom 2 throughout the 
process (see Fig. Pl. The total energy values along this 
unrelaxed lowest energy path are given in Fig. [3 (crosses), 
and the maxima of these curves give migration barriers 
Qb of 2.4 and 1.1 eV for the Ga -> VGa and Mnca ^ Vca 
jumps, respectively. The migration barriers for Mn/Ga 
migration in Fig.|21are qualitatively different, but similar 
barrier shapes with one and two saddle points have been 
observed for Ga 0] and substitutional Si iJ/S] migration 
in GaAs, respectively. Assuming that Qb is not affected 
by the presence of the Mnca atom, the activation bar- 
rier for pair dissociation may be estimated as Qt + Eb, 
yielding the dissociation barrier of 2.2 eV. 

We also considered another scenario within the 
formation-dissociation diffusion mechanism, where the 
Mn-VGa pair diffuses as an entity. The MnGa could ki- 
netically trap vacancies enabling the pair to move for 
some sufficiently long distance before dissociating, if the 
migration barrier for the vacancy to jump past the Mnca 
to another neighboring site were significantly lower than 
the MuGa^VGa pair dissociation barrier. However, for 
this jump we find a non-relaxed activation barrier of 
Qb — 2.3 eV, that is virtually the same value as for the 
pair Mnca-VGa dissociation. Thus such correlated diffu- 
sion will not occur and requires no special consideration. 

We improve the estimated migration barriers by al- 
lowing constrained relaxations as follows. The migrating 
Ga/Mn atom is allowed to relax along the [001] axis, the 
As atoms 1-3 in Fig. ^ are allowed to relax in the (110) 
plane, and the As atoms 4-7 are allowed to relax in the 
(110) plane. This calculation is carried out for the ini- 
tial and saddle point configurations of the non-relaxed 
calculations. For Mn we also calculate the local min- 
imum point at the center of the migration path. The 
resulting energies are given in Fig. [3 (circles), and we ob- 
tain the improved barrier maxima of 1.6 eV and 0.8 eV 
for the Ga and Mn migration, respectively. Although 
the Mn rich metallic environment studied in this work 
differs from pure GaAs, the migration barrier for the 
vacancy mediated Ga self-diffusion of 1.6 eV is well in 
agreement with the first principles calculation for a neu- 
tral vacancy |l8j] . and also with the experimental values 
around 1.5-1.9 eV ^ ISl IS^ in pure GaAs (Table lU. 
We wish to remark that even the improved migration 
barriers given above may be overestimated because only 



TABLE I: Calculated migration barriers Qt for Ga self- 
diffusion and Mn diffusion via the Ga vacancy mechanism. 
The and Q" are the values from calculations with no re- 
laxation and constrained relaxation, respectively. "PWPP" 
denotes the plane wave pseudopotential method, and "Expt" 
experiment. 



Qt' (eV) Ql' ( eV) PWPP (eV) Expt (eV) 



Process 

Ga^ Vca 2.4, 2.3" 1.6 
MUGa ^ Vca 1.1 0.8 
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° Ga self-diffusion around a Mnn,. ^ p^. |2?il2lll23|. 



nearest neighbor atoms are allowed to relax. Neverthe- 
less, considering the calculated barriers, the probability 
of the MnGa~VGa exchange vs. the Vca migration jump 
differs at relevant temperatures by a factor of 10^. For 
the dissociation of a MnGa~VGa we obtain an activation 
barrier of 1.4 eV by combining the VGa migration bar- 
rier with the MnGa~VGa negative binding energy. The 
Boltzmann factor for this dissociation barrier is 100 times 
larger than that for the VGa migration barrier, but still 
10^ times smaller than the MuGa-VGa exchange proba- 
bility. This shows that the bottleneck for MuGa diffusion 
via the pair formation-dissociation mechanism is the mo- 
bility of VGa- However, at large MnGa concentrations, 
MuGa^VGa pairs are formed more frequently, and due to 
the low pair dissociation and exchange barriers the VGa 
diffusivity is expected to increase. 

The calculated binding energies ^2§\ and the migration 
energies Q" (Tabled are used to study structural evo- 
lution with lattice kinetic Monte Carlo simulations using 
the Casino-LKMC code |2i|. (Ga,Mn)As is studied in 
a supercell of more than 100 000 Ga atoms, where ran- 
domly picked Ga sites are either replaced with Mn repre- 
senting a Mn concentration [Mn] of 5 % or 0.5 % (~ 10^^ 
or 10^° cm~'^), or left empty representing a Ga vacancy 
concentration of 10"'^* cm~^, mimicking the experimental 
samples of Ref. . Here [Mn] is defined as the propor- 
tion of the number of Mn atoms to the number of Ga 
sites. 

We study the clustering due to Mn redistribution by 
evaluating the number of Mn atoms in clusters as a 
function of annealing time and temperature, given in 
FigsEl(a) and (b). At [Mn] of 5% and 200°C [Fig.|Sl(b)] 
we find practically no clustering, which is fully consistent 
with the experimental finding that Tc of samples with 
[Mn] = 6% remains practically unchanged during long- 
time annealing {t > 10 h) at temperatures lower than 
190°C H. The Mn-Mn attraction combined with efli- 
cient Mn diffusion leads to increasing clustering with in- 
creasing temperatures from 250°C on [FigsOKa) and (b)]. 
The clustering rate exhibits a typical Arrhenius depen- 
dence, i.e. the logarithm of the clustering rate depends 
linearly on the inverse temperature. We find further that 
an increase in Mn concentration (at fixed temperature) 
increases the clustering rate as well. This increase is due 
to the fact that fewer migration jumps are required for 
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FIG. 3: Percentage of Mn atoms included in clusters as a 
function of annealing time for the Mn concentrations of 0.5 % 
(a) and 5 % (b), and cluster portions for the Mn concentration 
of 5 % as a function of time at temperatures of 250° C (c) and 
300° C (d). 



Mn to reach another Mn atom. At the same time signif- 
icant clustering starts to occur at lower and lower tem- 
peratures. This behavior is seen explicitly in Figs|31(a) 
and (b), where the clustering rate for [Mn] = 5% at 250 
(300)°C is approximately the same as for [Mn] = 0.5% at 
300 (350)°C. The clustering rate changes rapidly along 
the temperature, as seen for [Mn] = 5% in Figs El (c) 
and (d), where the temperature is increased from 250 to 
300°C. We get a similar increase in clustering rate by 
increasing the Mn concentration from 5% to 8% e.g. at 
250°C. The largest cluster size after annealing at 250°C 
for 24 h at [Mn] = 5 and 8% is 14 and 26 Mn atoms, 
respectively. The increasing number of large clusters 
may further indicate the formation of a secondary MnAs 
phase, as observed in growth of (Ga,Mn)As samples with 
Mn concentrations beyond 7 % 

The Curie temperature Tc depends mainly on the con- 
centration of Mn clusters [cl] [defined as the proportion 
of the number of clusters (including monomers) to the 
number of Ga sites] |llL Hq . i.e. Tc remains approxi- 
mately unchanged at fixed [cl] even when [Mn] is varied. 
Further, in the mean field approximation the concentra- 
tion dependence of Tc is linear, i.e. Tc cx[cl]. At the 



Mn concentration of [Mn] = 5% we find that [cl] evolves 
as follows: initially [cl] is 3.6%, and after annealing at 
250°C or 300°C for 24 h [cl] drops to 3.3% or 1.5%, re- 
spectively. In Ref. ^3 show that [cl] values of 6.3 and 
3.1 % correspond to Tc values of 660 and 220 K, respec- 
tively, and using the mean field approximation for Tc, 
the drop in [cl] from 3.6 to 3.3% at 250°C corresponds 
to a drop in Tc of 40 K. Experimentally, annealing of 
samples with [Mn] = 6-8% at temperatures from 215 to 
250°C for 24 h reduces Tc by 10-50 K 0,13, to which 
our result is in close agreement. Furthermore, Stanciu et 
al. find that annealing at 275°C induces a rapid drop in 
Tc of 20 K in 4 h 8] , which is bracketed by our estimated 
drops in the same time at 250°C and 300°C. 

To conclude, we have studied the redistribution of sub- 
stitutional Mnca in (Ga,Mn)As. The binding energies 
and migration barriers are calculated from first princi- 
ples. We show that substitutional Mn atoms use Ga va- 
cancies as a vessel for diffusion. Lattice kinetic Monte 
Carlo simulations yield that at annealing temperatures 
above 250°C the substitutional Mn redistribute into clus- 



ter configurations, and that the clustering rate increases 
along temperature and/or Mn concentration [Mn]. This 
clustering on the other hand lowers the Curie tempera- 
ture Tc, thus explaining quantitatively the drop in Tc 
observed in long-time annealing experiments. The in- 
creased clustering seems to impose a fundamental limit 
on Tc- Further, at large [Mn] beyond 8% clustering takes 
place already during sample growth, and formation of 
large Mn clusters indicates precipitation. 
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